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Abstract 


A  new  representation,  called  the  total  hazard  construction,  of  dependent 
random  variables  by  means  of  independent  exponential  random  variables  is 
introduced.  Conditions  which  imply  association  of  nonnegative  random 
variables  are  found  using  this  construction.  Furthermore,  new  conditions 
whichjmply  stochastic  ordering  between  two  nonnegative  random  vectors  are 
obtained.  These  strengthen  previous  results  of  the  authors.  Further 
applications  in  reliability  theory  and  in  simulation  are  indicated. 

AMS  subject  classification:  Primary:  60E05,  Secondary:  60K10,  65C05. 

Key  words  and  phrases:  hazard  functions,  exponential  distribution, 
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load-sharing,  multivariate  Pareto  distribution,  MIFRA  distributions,  variance 
reduction  in  simulation,  imperfect  repair. 


1.  Introduction 


Let  T  be  an  absolutely  continuous  random  variable  to  be  thought  of  as  a 
lifetime  of  a  device.  If  F  and  f  are,  respectively,  the  distribution  and 
the  density  functions  of  T  then  7=1-  F,  a  =  -  log  7  and  x  =  f/7  are, 
respectively,  the  survival,  hazard  and  hazard  rate  functions  of  T.  It  is 
easy  to  verify  that  if  X  is  standard  (that  is,  mean  1)  exponential  random 
variable  then 

(1.1)  T  =  A_1 (X ) 

sati sfies 

(1.2)  f  ^  T 

where  ,-tl  denotes  equality  in  law. 

Equations  (1.1)  and  (1.2)  suggest  that  study  of  A  can  shed  light  and 
understanding  on  the  stochastic  behaviour  of  T.  One  purpose  of  this  paper  is 
to  introduce  and  study  a  multivariate  analog  of  (1.1)  and  (1.2).  This  is  done 
in  Sections  2  and  3  where  it  is  shown  how  to  transform  a  set  of  independent 
standard  exponential  random  variables  Xj,...,X  into  a  random  vector  _T 
which  is  stochastical ly  equal  to  a  given  random  vector  T.  The 
transf ormat i on  uses  multivariate  hazard  functions  (to  be  defined  in  Sections  2 
and  3)  and  will  be  called  the  total  hazard  construction. 

In  Section  4  we  list  conditions  under  which  the  transfornati on 
(Xj,...,X  )  -*■  T  is  increasing  (in  this  paper  'increasing'  stands  for 
'nondecreasing'  and  'decreasing'  for  'nonincreasing' )  in  each  of  the  X^'s. 

A  A 

Under  these  conditions,  then,  the  random  variables  Tp...,7n  (and  hence 


2 


Tj,...,Tn)  are  associated  in  the  sense  of  Esary,  Proschan  and  Walkup  (1967). 

Let  S  and  T  be  two  absolutely  continuous  random  variables  with  hazard 
functions  Q(w)  =  -  log  P{S>w}  and  R(w)  =  -  log  P{T>w}.  Then,  using  the  same 
standard  exponential  random  variable  X,  one  can  apply  (1.1)  and  (1.2)  to 

A  A 

obtain  S  and  T,  defined  on  the  same  probability  space,  such  that 

(1.3)  S  =  Q-1(X)  S,  T  =  R"1 (X)  T. 

If 

(1.4)  Q(w)  >  R(w),  w  >  0, 

r  f  A  A  c  t 

then,  from  (1.3),  S  -  S  <  T  -  T.  Thus  we  see  that  (1.4)  implies  stochastic 
ordering  of  S  and  T.  In  Section  5  we  obtain  a  multivariate  extension  of 
this  result,  again  using  the  total  hazard  construction  of  Sections  2  and  3. 
Further  applications  are  given  in  Section  6. 

A  random  variable  S  is  said  to  be  stochastically  smaller  than  a  random 
variable  T  (denoted  S  T)  if  P{S  >  u}  <  P{T  >  u)  for  every  u.  A 
random  vector  SI.  =  (Sj,...,S  )  is  said  to  be  stochastically  smaller  than  a 
random  vector  T  =  (Tj,...,T  )  [denoted  S^  JT]  if 

(1.5)  g(S)  ^  g(T) 

for  every  increasing  Borel  measurable  real  function  g.  A  function  g  is 
called  increasing  if  (xj,...,xp)  ^  (yj,....yn)  implies 
g(x1,...,xn)  s  g(y1,...,yn)  where  (xlt...,xn)  <  (y x . ,yp )  means 
x^  <  y^,  i  =  l,...,n  .  It  is  well  known  that  S.  ^  T_  1  f  ancl  only  if 


3 


Eg(S)  <  Eg(T) 

for  every  increasing  Borel  measurable  real  function  g  for  which  the 
expectations  exist.  Also  S_  X  if  and  only  if 

P{S  eU)  <  p{T  e  U} 


for  every  Borel  set  U  which  has  an  increasing  indicator  function. 

2.  The  total  hazard  construction:  bivariate  case. 

Consider  two  nonnegative  random  variables  Tj  and  T£  with  absolutely 
continuous  joint  distribution  function  F,  joint  density  function  f  and 
joint  survival  function  7  defined  by  F(tpt2)  =  P{Tj  >  tp  T,,  >  t2) .  The 
conditional  hazard  rate  of  Ti  at  time  t,  given  that  T3_^  >  t  is  defined 
as 


(2.1)  A .  ( t )  =  lim  -j-  P{t  <  T.<  t  +  At|T  >  t,  T?  >  t} 
At+0 


tj-tj-t)'  1  *  °>  1  ’  1,2‘ 


Given  that  T2  =  t£»  the  conditional  hazard  rate  of  T^  at  time  t  >  t2  is 
defined  as 


(2.2)  X1(t|t2)  =  lim  P(t  <  Tj  <  t  +  At|Tj  >t,  T2  =  tg} 


At+0 


f  (t,t2) 


,  t  >  t 


2  * 


I 


(Substituting  t  and  t 2  in  the  right  hand  side  of  (2.2)  may  yield  0/0. 
Here  and  in  the  remainder  of  the  paper  such  a  ratio  is  interpreted  as  0. ) 
Similarly  define 


,  f(ti,t) 

(2.3).  X2(t|tj)  5 


Stj  F(tl,t) 


»  t  >  t.  • 


The  total  hazard  accumulated  by  time  t  by  the  random  variable 
given  that  minfTj^)  >  t,  is  defined  by 


(2.4)  Ai  (t )  =  / q  Xi(u)du,  t  >  0,  i  =  1,2. 


Given  that  T2  =  t2  and  that  Tj  >  t2,  the  total  hazard  accumulated  by  the 
random  variable  Tj  during  the  time  interval  [t2,t2+t)  is  defined  by 


t2+t 


(2.5)  Aj(t|t2)  =  /t  X1(u|t2)du,  t2  >  0,  t  >  0. 


Similarly 


t,+t 


(2.6)  A2(t|tj)  =  /t|  A2(u|t1)du,  tj  >  0,  t  >  0. 


Note  that  with  this  definition. 


(2.7)  -~r  Ai  (t|  t2  )  =  (t 2+t|  t2),  t2  >  0,  t  >  0, 


(2.8)  —  A2  (t )  1 1 )  =  X2(t1+tjt1),  tj  >  0,  t  > 


0. 


The  total  hazard  accumulated  hy  T^  by  the  time  it  failed  is  defined  as 


tkjy'Iw 


Ai (Ti )  if  <  T 2  and  as  A^(T2)  +  (T^ -T2 1  T2 )  ^  >  ^2*  Similarly 

define  the  total  hazard  accumulated  by  Tg  by  the  time  it  failed. 

Define  the  inverse  functions 

A^(x)  =  inf{t  >  0:  A  -  (t )  >  x} ,  i  =  1,2,  x  >  0, 

A^’1(xjt2)  =  i nf { t  >  0:  A1(tjt2)  >  x},  x  >  0,  t2  >  0, 

A21(xjt^)  =  i  n  f  { t  >  0:  a2  (t  1 1  ^ )  >  x},  x  >  0,  t1  >  0, 

2  2 

and  consider  the  functions  a^:R+  R  and  a2:R+  -*•  R  define  as  follows:  On 
2-1  -1 

{(Xj,x2)  e  R+:  A^  (xj)  <  A-?  ( x 2 ) } » 

(2.9)  a1(x1,x2)  =  Aj^Xj), 

(2.10)  a2(x ^,x2)  =  A^  (x^)  +  A2  (x2~A2(A^  (x^))|a^  (x^)), 

2  -1  -1 

and  on  {(Xj.Xg)  e  R+:  A^  (x^)  >  A2  (x2)} 

(2.11)  a1(x1,x2)  =  A2X(x2)  +  a;1(x1-A1(A21(x2))|a21(x2)), 

(2.12)  a^(xj,x2)  =  A2  ( x 2 )  . 


Motivated  by  the  fact  (see  Remark  2.2  below)  that  the  total  hazards 
accumulated  by  Tj  and  T£  by  the  time  they  failed,  are  independent  standard 
exponential  random  variables  we  introduce  the  following  total  hazard 
constructi on : 

Let  Xj  and  X2  be  independent  standard  exponential  random  variables 
and  consider  the  following  transformation: 


Theorem  2.1.  Let  (Tj.Tg)  be  defined  as  in  (2.13)  where  Xj  and  X£  are 
independent  standard  exponential  random  variables.  Then 

(2.14)  (Tltf2)  =*  (Tj,T2)  . 

Proof.  Let  A^(*|t2),  and  so  on,  be  the  analogs  of  a..,  Aj(»|t2),  and 

a 

so  on,  defined  in  (2.1)  -  (2.3)  with  T.  replacing  Ti#  i  =  1,2.  We  will 
show  that 

(2.15)  A^t)  =  Aj (t ) ,  t  >  0, 

(2.16)  A2(t)  =  A2(t),  t  >  0, 

(2.17)  A1(t|t2)  =  A 1  ( 1 1 12  ) ,  t  >  t2  >  0, 

(2.18)  A2(t|t1)  =  A  2  ( t  j t^  ),  t  >  t^  >  0, 

and  the  result  then  follows  from  the  fact  that 
A  x  ( • ) ,  A2(0,  A :  ( •  1 1 2 ) ,  A2(.|t1),  tj  >  0,  t2  >  0  uniquely  determine  f 
(see,  for  example,  the  explicit  formulas  in  Cox  (1972)  or  in  Shaked  and 
Shanthikumar  (1984a))  and  that 

AAA  A 

A^  (• ),  A2(.),  A1(*jt2),  A2|.(t1),  t1  >  0,  t2  >  0,  determine  the  joint 

A  A 

density  of  T^  and  T2  in  a  similar  manner. 

To  show  (2.15)  notice  that 

(fj  >  t,  t2  >  t  }  «  { aJ1(x1 )  >  t,  a‘1(x2)  >  t}  =  (x1  >  \l(t),  x2  >  A 2 ( t ) } 
Hence,  for  t  >  0, 

P(t  <  Tj  <  t  +At|f1  >  t,  T2  >  t} 

=  P{t  <  A^1 (X1 )  <  t  +  At|X1  >  AT(t),  X2  >  A 2 ( t ) } 

=  P { A  x ( t )  <  X1  <  Aj(t  +  At )  j  X  L  >  A 1  ( t ) ,  X2  >  A  2 ( t ) } 


=  P{X1  <  Aj(t  +  At)  -  A  j ( t ) } 

=  1  -  exp{-(A1(t  +  At)  -  Aj^(t))}. 


Thus 


X 2 (t )  =  Tim  P{t  C  Tj  <  t  +  At|fj  >  t,  Tg  >  t} 

At  +  O 


.  1 


lim—  (l-exp{-(A,  (t+At)  -  A-i(t))}) 
At+0A1: 


dt  Aj(t)  =  Aj(t) 


and  (2.15)  follows.  The  proof  of  (2.16)  is  the  same. 
To  prove  (2.17)  note  that  for  t  >  tg  >  0, 


{^1  >  t>  ^2  =  t2 1 

=  {A^1  (X2)  +  A‘1(XrA1(A21(X2)))A21(X2))  >  t,  A^1  (X2)  =  tg} 
=  (t2  +  A^Xj-AjUgJltg)  >  t,  A^(X2)  =  t2} 

=  { X ^  ^  A ( 1 2 )  +  A 2 ( t ~t 2I  tg),  Ag  (X 2 )  =  1 2 }  • 

Thus,  for  t  >  t2  >  0, 

Pft  <  T1  <  t  +  At J  T j  >  t,  fg  =  tg} 

=  P{t  <  A^(Xg)  +  A;1(X1-A1(Ag1(Xg))|  Ag1(X2))  <  t  +  At 
|  X 1  >  A^(tg)  +  A 2  ( t - 1 2 j  tg) ,  A 2  (X^)  =  1 2} 

=  P{  t-tg  <  A2  ^  ( X2  —A 2  ( 1 2  )|  tg)  <  t  +  At  -  tg 

|Xj  >  Ax(tg)  +  A  2  ( t  —  t  g  j  tg),  A  g 1  ( X  g )  =  tg} 

=  P  { Ax  (t  2)  +  A  j  (t-tg|  t  g)  <  X1  <  A  j  (t  2)  +  A  2  (t -t  g+At  |  tg) 
j  X 1  >  A^(tg)  +  A  2  ( t  — t  g|  1 2 ) »  Ag  (Xg)  =  tg} 

=  P  {  X  2  <  A 1  (t-t  g+At  |  t  g )  -  A  j  ( t  —  t  g  |  tg)}  . 


Hence,  for  t  >  >  0, 

X,(t|t2)  =  lim  ^r{l-exp{-(A1(t-t2+At|t2)  -  Aj  (t-t2|  t?) )} ) 

At+O 

=  ft  Ai  (t"t2l  t2)(2=?)  xl(tlt2) 

and  (2.17)  follows.  The  proof  of  (2.18)  is  the  same,  n 

Remark  2.2.  From  (2.13)  one  easily  obtains:  On  { T ^  <  T2|  ,  =  a^Tj)  , 

X2  =  ^2(T^ )  A2 ^2~^1 1  ^1  snd  on  {T^  >  T2)  , 

X1  =  Ai  (Tp )  +  (Tj-T2|  T2 )  ,  X2  =  a2(T2)  .  Thus  Xj  and  X2  are  the  total 

A  A 

hazards  accumulated  by  Tj  and  T2  by  the  times  they  failed.  Since  (Tj, 

T2)  -t  (Tp  T  ),  it  follows  that  the  total  hazards  accumulated  by  T1  and 

T2  by  the  time  they  failed  are  independent  standard  exponential  random 

variables  (see  also  Remark  3.2). 

Examp1 ~  1  (Freund  (1961)  distribution).  The  Freund  distribution  provides  a 
model  in  which  the  exponential  residual  life  o *  one  component  depends  on  the 
working  status  of  another  component.  It  has  the  density 

(2.19)  f  (t  1 ,  t2  )  =  a8  'exp{  -(a+3-0  '  )t  ^  -  il't^  if  0  <  tj  <  t?  , 

=  a 1 8  expt-a'tj  -  (a+8-a')t2)  if  0  <  t/>  <  , 


where  a,  a',  3,  3 


are  nonnegative  parameters.  For  this  distribution 


A 1  ( t j  1 2 )  =  ot't,  t  >  0,  t2  >  0  , 
A2(tjt1)  =  g't,  t  >  0,  tj  >  0  . 


If  (T ^ , T2 )  has  the  joint  distribution  (2.19)  then  from  (2.13)  and  (2.14) 


it  follows  that 

it  has 

the 

same  distribution 

as 

(t1,t2) 

where 

(2.20)  fx  = 

a  ^X  ^ 

if 

a  *Xj 

<  6 -1X  2  , 

= 

6_1X2  + 

(a') 

■1(X1-8"1aX2) 

if 

-lv 

a 

C\J 

X 

H 

1 

CD. 

A 

f2  = 

a  ^X^  + 

on 

■1(X2-ct"16X1) 

if 

a  1X1  < 

;  e_1x2  , 

6_1X2 

if 

a  1X1 

>  8  _1X  2  , 

and  X1  and  X2  are  independent  standard  exponential  random  variables. 
Representati on  (2.20)  can  be  rewritten  as 


(2.21) 

T1 

-  "ly 

-a  a  ^ 

if 

a 

-ix 

xi 

< 

b'1x2  , 

=  (a  '  )  ^Xj 

+  B_  1  ( 1 

1 

X 

ro 

if 

a 

-ix 

xi 

> 

b_1x2  , 

f2 

CM 
X 
» — w 

1 

CQ 

11 

+  a_1(l 

-  Kn, 

IT'  1 

if 

a 

-ly 

xi 

b_1x2  . 

=  6_1X2 

if 

ct 

-ix 

xi 

> 

b_1x2  . 

Representation  (2.21)  is  identical  to  (7)  of  Shaked  (1984). 

The  example  will  be  continued  later. 

Example  2  (bivariate  Pareto).  Let  (T1,T2)  have  the  joint  survival  function 

F(t1,t2)  =  (1+tj+tj,)-1  ,  tj  >  0,  t2  >  0  . 


It  is  not  hard  to  verify  that  in  this  case 


(2.22)  A^(t)  =  A2 (t: )  -j  log  (l+2t),  t  >  0  , 

(2.23)  A1(t| t2)  =  2  log  (1+  j^),  t  >  0,  t2  >  0  , 

(2.24)  A2(t|  t x)  =  2  log  (1  +  j^),  t  >  0,  tj  >  0  . 

Some  algebra  shows  (using  (2.13)  and  (2.14))  that  (Tj.Tg)  has  the  same 
distribution  as  (T^.T^)  where 

«  1 

T1  =  7  (e  -1)  if  Xx  <  X2  , 

2X 

=  exp{^-X1  +|x2}  -  |(e  2+l)  if  ^  >  X?  , 

2  X 

=  exp {?X ^  +  "2^2^  “  Tp  +1)  if  ^1  *  ^2  ’ 

1  ^X« 

=  -g  (e  -1 )  i  f  X^  >  X^  ; 

here  X^  and  ^  are  independent  standard  exponential  random  variables. 

This  example  will  be  continued  later. 

3.  The  total  hazard  construction:  multivariate  case. 

Consider  a  random  vector  T^  =  (T^,...,T  ),  n  >  2,  with  absolutely 

continuous  joint  distribution  function.  In  this  section  we  describe  the  total 

hazard  construction  of  a  random  vector  T  =  (Tj,...,T  )  such  that  1  ^  T  , 

The  construction  will  be  described  in  n  steps  numbered  1  through  n.  In  Step 

1  an  index  ji  is  chosen  at  random  from  {1,2, ...,n)  and  then  T.  is 

J1 

determined.  Upon  entering  Step  k,  2  <;  k  <  n,  the  random  variables 

A  A 

T.  ,...,T.  have  already  been  determined  where 

J1  Jk-1 

J  =  C  { 1, . .. ,n}  .  In  Step  k  an  index  jk  is  chosen  at 


random  from  J  =  -  J  and  then  T,  is  determined. 

Jk 

We  need  to  extend  and  slightly  modify  the  notation  of  Section  2. 

For  J  =  C  {l,...,n}  let  tj  denote  (tj  ,...,tj  ).  If 

7  =  then  tj  denotes  (ti^,...,ti  ^_).  Let  e  =  (1,...,1). 

The  length  of  e  will  vary  from  one  formula  to  another,  but  it  will  be  always 

possi.ble  to  determine  it  from  the  expression  in  which  e  appears. 

For  J  C  {1,  ...,n}  and  i  e  U  let  x^  (tj  Tj  =  tj,  T-j  >  te)  denote  the 

conditional  hazard  rate  of  at  time  t  given  that  =  tj  and  that 

T-t  >  te  where  t  >  V  t.  =  max{t.:j  e  J}  .  If  J  =  0  then 

— —  ^  i  J  J 

jeJ 

V  t.  =  0  .  Formally,  for  i  e  IT  , 
jeJ  J 


(3.1)  ^-j  (^ | Xj  =  Ij  >  te) 

=  lim  P{t  <  T.  <  t  +  Atj  T . 

AUO  At  1  '-0 


tj,  T-r  >  te),  t  >  V  t. 
^  ^  "  jeJ  J 


(J  may  be  empty).  The  absolute  continuity  of  T_  ensures  that  this  limit 

exists.  To  save  space  we  sometimes  supress  the  condition  Tj  >  te  and  just 

write  x  ^ ( t  j  Ij  =  but  the  reader  should  keep  in  mind  that  '•  '  means 

T-y  >  te  with  t  being  the  same  as  the  first  argument  of  x.  .  The  function 
— 0  —  i 


Xi (.j Tj  =  tj,  •)  will  be  of  interest  for  us  only  on  the  (random)  interval 

(max  T.,  min^T.],  however,  to  avoid  a  discussion  of  such  random  hazard  rate 
j  e  J  J  i  e  J 

functions  [such  a  discussion  can  be  found  in  Arjas  (1981)]  we  do  not  emphasize 

this  point  here.  Note  however  that  x^(t|Tj  =  t^,  •)  is  well  defined  for 

every  t  >  V  t,  . 

jeJ  J 

For  i  e  7  the  total  hazard  accumulated  by  T^  during  the  time  interval 

[V  t.,  V  t.  +  t),  t  >  0,  is  defined  by 
jeJ  J  jeJ  J 

V  t  .+t 
jeJ  J 

A*  (t|  T ,  =  t,)  =  /  A  •  (u|  T  i  =  t,,  •  )du ,  t  >  0,  i  e  J  . 


(3.2) 
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When  J  =  0,  A^(tjTj  =  tj)  will  be  simply  denoted  by  A^(t)  . 

We  will  introduce  now  a  notation  for  the  total  hazard  accumulated  by 

Ti  by  time  t.  Fix  t  >  0  and  suppose  that  it  is  given  that  T. 

'  J1  J  k -1 

(k  >  1)  failed  at  times  t^  .....t^  ,  respectively  (t.  <...<t.  <  t) 

J1  Jk-1  Jk_1 

and  that  all  the  other  Tj's  are  alive  at  time  t.  For  i  ^{jj,...,jk  ^1, 


denote 


hli  -j  (t|t-  ,...,t-  )  =  A-(t.  ) 

M  J  i . Jk-i  J1  Jk-1  1  J1 

k-i 

+  I  A  •  (t .  t  ■:  |  T  <  =  t .....  ,T  •  t .  ) 

1=2  1  h  J n-l|  J1  J1  jji-1  J£-l 

+  A.j(t-t.  T.  t.  t.  )  . 

1  Jk-1  J1  J1  Jk-1  Jk-1 


Also  denote  (corresponding  to  the  case  k  =  1) 


(t )  =  Ai  (t),  t  >  0. 


Note  that  for  tj  <  tj  and  i  j;  { 3V . . . ,  jk_1> 

^i|  •;  ■:  -j  1 1  ,-♦•••  »t  •  )  “  A .  (t  .  ) 

1 1 J 1 , - - . , Jk_l  Jk_l  Jl  Jk_l  1  Jl 


+  1  t  •  |  T  ,•  t T .  t  •  )  . 

1=2  JZ  Jn-1  J1  J£-l  Ji-1 


The  total  hazard  accumulated  by  T-j  by  the  tine  it  failed,  given  that 

Ti  was  the  k-th  Tj  to  fail  and  that  Tj^,...,Tj  failed  before  T^ ,  is 

.  .  (T,|T  .  T  )  . 

^  J  i » •  •  •  *  ^k-1 
Define  the  inverse  functions 


*  *  *  “  ***  •*  '**  * »*  *•*  *»*  *«“ 


x  >  0, 


A^(x)  =  i  nf  { t  >  0:A.(t)  >  x) ,  i  =  1,2,. ...n, 
and  for  nonempty  J  C  (1, . . . ,n} ,  tj  >  0e  and  i  t  J, 

Ai 1  ^x|lj  =  tj  )  =  i n f { t  >  0:  a.  (t|Tj  =  )  >  x},  x  >  0. 

Motivated  by  the  fact  (see  Remark  3.2  below)  that  the  total  hazards 
accumulated  by  Tj,...,Tn  by  the  time  they  failed,  are  independent  standard 
exponential  random  variables,  we  introduce  and  study  the  following  total 
hazard  construction: 

Let  X1,...,Xn  be  independent  standard  exponential  random  variables. 


Step  1.  Let  j]_  be  the  (random)  index  (which,  by  absolute  continuity,  is 

unique  with  probability  1)  such  that  aT1  (X.  )  =  minuj1  (X . ):  1  =  1,  ...,n} 

Jj  Jj  >  ' 

and  define 


(3.3) 


h 


=  aTX(x  ) 
J1  J1 


Step  k  (k  =  2,...,n).  Given  that  Steps  l,2,...,k-l  resulted  in 

A  A 

T-j  =  *  •  •  •  »Ti  =  t-  let  J  =  { ,. .. ,  ju  i) .  Let  ju  be  the 

J1  J1  Jk-1  Jk-1  * 

(random)  index  (which,  by  absolute  continuity,  is  unique  with  probability  1) 

such  that 


^  i  ♦  -j  I  -j  i  ■:  1 1  •,...,  t .  )  |  T  .=  t T  .  =  t  . 

Jk  Jk  Jki  Jl’ •**’Jk-l  Jk-1  J1  Jk-1  J1  J1  Jk-1  Jk-1 

mi  n  { A^  [X^  ^-jl  i  ;  (f .  1^,-  )  |  ^  .■  L  ■  ,...,T-  =  t- 

ie7  M J  i» - Jk-i  Jk-1  J 1  Jk-1  J 1  J1  Jk-1  Jk- 


-N  V  v~  C  s' 


1  •'wV.'i  •.*  O 0  V- 
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It  is  easy  to  verify,  by  induction,  that  the  arguments  of  a7*  and  a7*  ,  in 

1  Jk 

the  above  expression,  are  nonnegative.  Having  chosen  the  (random)  index 
as  described  above,  define 


(3.4.1)  T.  =  T. 

i  Jk  Jk-1 

+  A/  [X.  -  h  i  (Ti  |T.  ,...,Ti  )|Ti  «  t  -  t  ] 

Jk  Jk  Jkl Jl’***’Jk-l  Jk-1  Jk-1  J1  Jk-1  Jk-1 


More  expl  icitly. 


(3.4.1  i )  T.  =  T  +  aTV  -  a.  (T  ) 

Jk  Jk-1  Jk  Jk  Jk  J1 

k-1 

-I  A,  (fj  -f,  |Ti  =  T.  =T  )|  T  =  T . T  =T  ]. 

1=2  Jk  Jz  jji-1  J1  J1  Ji-1  J£-l  J1  J1  Jk-1  Jk-1 

For  example,  if  n  =  3  and  Aj1^)  <  A^(X 2 ),  A j1  (X x)  <  Aj1^)  an< 
[X2"A2(A^  (Xj))|T^  =  A^  (X^)]  <  A3  [X^  ”  A^(a^  (X1))|T1  =  A^  (X^)]  then 


(3.5)  Tj  =  Aj  (Xj)  , 

(3.6)  T2  -  fj  +  A21[X2  -  A2(T1)|T1  =  T1] 

+  a2  fX2  '  *2|l(fl!V!Tl  =  ’ 

(3.7)  T3  =  T2  *  A3V3  -  A^)  -  a3^2  -  fllTl  =  =  fP  T2 

A  _  J  AAA  *  A 

T2  +  A3  fX3  “  *3 | l,2^T2l T1,T2 H  T1  =  Tl»  T2  =  T2^  * 


Theorem  3.1.  Let  _T  be  as  defined  in  (3.3)  and  (3.4)  where  Xj,...,X  are 
independent  standard  exponential  random  variables.  Then 


(3.8)  T  -  T  . 


To  prove  Theorem  3.1  define  analogously  to  the  definition  of  in 

A 

(3.1)  with  J  replacing  Then  all  that  is  needed  to  complete  the  proof 

(see,  e.g.,  Shaked  and  Shanthikumar  (1984a))  is  to  show  that  for  all 

J  C  {l,...,n}  ,  i  t  J,  t  >  V  t.,  we  have 

jeJ  J 

(3.9)  X1  (t|  Tj  =  tj,  •)  =Xi(t|Tj  =  tj,  •)  . 

The  proof  of  (3.9)  is  similar  to  the  proofs  of  (2.15)  -  (2.18)  but  is 
notationally  more  involved.  We  omit  the  details. 


Remark  3.2.  It  can  be  shown,  using  (3.3)  and  (3.4)  that,  for  every 

A  A  A 

permutation  (j1,...,jn)  of  (l,...,n),  on  (T.  <  T.  <...<T.  }  we  have 

J1  J2  Jn 


=  *4  (T,  )  , 

J1  J1 


♦Hi  t  (Ti  lTi 
Jkl Jl >  —  »Jk-l  Jk  J1 


T.  ),  k  -  2, . . .  ,n  . 
’  Jk-1 


thus  Xj,...,Xn  are  the  total  hazards  accumulated  by  T^,...,T  by  the  times 

c  t  A 

they  failed.  Since  1_  -  ]_  it  follows  that  the  total  hazards  accumulated  by 
Ti»...,Tn  by  the  times  they  failed  are  independent  standard  exponential 
random  variables.  This  fact  generalizes  Theorem  2.2  of  Schechner  (1984).  It 
also  follows  from  Section  4.5  of  Aalen  and  Hoem  (1978)  and  Proposition 
2.2.11  of  Jacobsen  (1982). 


Remark  3.3.  It  should  be  emphasized  that  the  total  hazard  construction  (3.3) 
and  (3.4)  is  theoretically  and  practically  different  than  the  following  well 
known  standard  construction  (see,  e.g..  Law  and  Kelton  (1982),  p.  268  or 


Let  Uj,  U2,...,Un  be  independent  uniform  [0,11  random  variables  and  let 
X  =  )  be  an  absolutely  continuous  random  vector.  Define 

(3.10)  TJ  =  inf{tj :  P{T1  >  t^  >  l^}  , 

(3.11)  TjJ  =  inf {tk :  P{Tk  >  tk|  Tx  =  T{,...,Tk_1  =  Tk-1>  >  Uk) ,  k  =  2,3 . n. 

Then 

(3.12)  (T{,...,r)  &  (T^...,^). 

Although  the  construction  defined  by  (3.3)  and  (3.4)  is  different  than 
the  one  defined  by  (3.10)  and  (3.11),  the  results  which  follow  from  (3.3)  and 
(3.4)  have  analogs  which  follow  from  (3.10)  and  (3.11).  These  analogs  will  be 
noted  throughout  the  sequel. 

4.  An  application:  association  of  random  variables. 

4. 1.  The  bivariate  case. 

Let  T^  and  T2  be  nonnegative  absolutely  continuous  random  variables 

a 

as  in  Section  2  and  let  T^  and  T ^  be  defined  as  in  (2.13).  Since 
(Tp  T^)  (Tp  Tg)  we  will  not  distinguish  in  this  section  between  J  and 

A 

X  and  just  write  X  • 

If  the  functions  a^  and  a2  defined  in  (2.9)  -  (2.12)  are  increasing 

in  each  argument  when  the  other  argument  is  held  fixed  then  from  (2.13)  it 

follows  that  T^  and  T2  are  associated  in  the  sense  of  Esary,  Proschan  and 

Walkup  (1967).  Association  is  a  property  which  yields  important  probability 

inequalities  and  is  particularly  useful  in  reliability  theory  (see,  e.g., 


Barlow  and  Proschan  (1975)).  Thus  it  is  of  interest  to  find  conditions  which 


vvw 
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imply  that  and  &2  are  increasing. 

Theorem  4.1.  If  for  all  y  >  0,  u  >  0,  b  >  0,  i  =  1,2, 

(4.1)  A^y+u)  -  Ai  (y)  +  A^bjy+u)  <  A^  (u+b| y ) 
then  and  T2  are  associated. 

Remark  4.2.  Intuitively,  for  i  =  1  say.  Condition  (4.1)  says  that  the 

larger  T2  is  (compare  T2  =  y  to  T2  =  y  +  u)  the  smaller  is  the  potential 

hazard  that  can  be  accumulated  by  by  the  time  y  +  u  +  b  (see  Figure 
4.1).  Thus,  roughly  speaking,  the  larger  T2  is  the  larger  Tj  is  and  so 
the  association  of  Tj  and  T2  is  not  surprising.  For  a  similar  result  see 
Arjas  and  Norros  (1984). 

Proof  of  Theorem  4. 1.  We  just  have  to  show  that  ai(xi,x2)  and  a2(xj,x2) 
increase  in  x^  and  x2.  The  result  then  follows  from  Esary,  Proschan  and 
Walkup  (1967). 

Consider  a^(x^,x2).  It  is  easy  to  see  that,  for  a  fixed  x2,  a1(x1,x2) 

increases  in  x-y.  Thus  one  just  has  to  show  that,  for  a  fixed  xj,  the 

function  a^(x^,x2)  =A2*(x2)  +  a^ [x^-A^ (a2^ (x2 ))) a2^ (x2)]  increases  in 
x2  e  [0,  A2(A^1(x1))]  .  Denote  t  =  A^1^),  =  A^1  (x  1 )  .  We  need  to  show 

then  that,  for  each  ty  >  0, 


(4.2)  a(tj,t)  -  t  +  Aj1[ Ax (tx  )  -  Aj(t)|t] 


i ncreases  in  t  e  [0,  t, ]  . 


Rewrite  (4.2)  as  A^tj)  -  Aj(t)  =  A^aft^t)  -  t|t].  Similarly,  if 
u  >  0  and  t  +  u  <  tp  then 

Ai  (t  i )  -  A^t+u)  =  A1[a(t1,t+u)  -  t  -  ujt  +  U].  Thus 

(4.3)  Aj(t+u)  -  A^(t)  +  Aj[a(tpt+u)  -  t  -  u|  t  +  u]  =  A1[a(t1,t)  -  1 1 1 ] . 

But,  by  (4.1)  [here  LHS  =  left  hand  side,  RHS  =  right  hand  side], 

LHS  (4.3)  <  A1[a(t1,t+u)  -  t j  t ] .  Thus, 

A,[a(t1,t)  -  t|t]  <  A1[a(t1,t+u)  -  t j t ] .  But,  for  each  t,  Aj(b|t)  is 
increasing  in  b.  Hence  a(tpt)  <  a(tj,t+u),  that  is,  a(tpt)  increases 
in  t  e  [0,tj]. 

Similarly  it  can  be  shown  that  a2(xl»x2)  increases  in 
Xj  g  [0»  (A2  (X2 ))  ]  •  ii 

Condition  (4.1)  can  be  written  by  means  of  derivatives  (if  they  are  well 
defined ) : 

Theorem  4.3.  If  for  i  =  1,2, 

(4.4)  A.  (y)  +  |_.  a.  (bj  a)|  a=y  <  A-  (y+bjy ),  b  >  0,  y  >  0 

(provided  the  derivatives  in  (4.4)  are  well  defined)  then  and  T2  are 

associated. 

Proof.  Rewrite  (4.1)  as:  for  i  =  1,  2,  y  >  0,  b  >  0, 

(4.5)  [A.  (y+u)  -  A.  (y )]  +  [A.(bjy+u)  -  Ai(b|y)]  <  A^b+ujy)  -  Ai(b|y),  u  >  0. 


Dividing  (4.5)  by  u  >  0  and  letting  u  4-  0  one  obtains  (4.4)  from  (4.5).  To 


\  ^  i— y  rrrrF^V^TTr^ywro^ »»»»» rg^jnjrr. 


obtain  (4.5)  fi  "m  (4.4)  integrate  (4.4)  with  respect  to  the  dummy  variable 

y.  « 

Example  1  (continued).  From  (2.21)  it  is  easily  seen  that  if 

(4.6>  a  <  a'  and  B  <  81 

A  A 

then  and  T2  are  increasing  functions  of  X}  and  X2.  Hence  if 

(Tj.To)  has  the  Freund  distribution  with  parameters  satisfying  (4.6)  then 
Tj  and  T2  are  associated.  This  result  has  been  obtained  also  in  Shaked 
(1984). 

Example  2  (continued).  Differentiating  (2.22)  -  (2.24)  one  obtains 

Ai  (y )  =  i+2y  ’  ^  2’ 

Ai  ^+b|y)  =  l+2y+b  ’  1  1’2’ 

fa  Ai  (bla)la=y  =  ( l+2y )  ( l  +  2y+b)  ’  1  =  1,2‘ 

It  is  not  hard  now  to  verify  (4.4).  Hence  if  (T^,T2)  has  the  bivariate  Pareto 
distribution  then  Tj  and  Tg  are  associated. 

This  result  is  not  surprising.  Shaked  (1977)  has  shown  that  the 
multivariate  logistic  distribution  of  Malik  and  Abraham  (1973)  has  some 
positive  dependence  properties.  Since  the  multivariate  Pareto  distribution  is 
a  simple  transformation  of  the  multivariate  logistic  distribution  it  follows 
that  also  the  multivariate  Pareto  distribution  has  some  positive  dependence 
properties.  It  is  not  hard  to  find  other  representations  of  T^  and  J?  as 
incresing  functions  of  independent  random  variables. 


Example  3  (bivariate  Gumbel  exponential  distribution).  Let  Tj  and  T2  have 
the  bivariate  survival  function 

F(t1,t2)  =  exp{-t1-t2-0t1t2),  t1  >  0,  t2  >  0, 

where  9  e  [0,1]  is  a  fixed  parameter.  Here 

(4.7)  Aj(t)  =  A2(t)  =  t  t2  ,  t  >  0 

(4.8)  Aj  (u|  t )  =  A2(u|t)  =  u  +  Qtu  -  log  ,  t  >  0,  u  >  0. 

The  inverses  Aj1  and  A^  do  not  have  as  simple  expressions  as  in  Examples 
1  and  2  but  it  is  still  not  hard  to  check  (4.4).  Differentiating  (4.7)  and 
(4.8)  one  obtains 

Ai  (y)  *  1  +  ey  ,  1 
X-  (y+b|y)  =  1  +  ey  - 
h  A,(b|a)|a,y  =  eb 

Substituting  these  in  (4.4)  it  is  seen  that  (4.4)  does  not  hold.  Thus  we 
cannot  show  that  Tj  and  T£  are  associated.  In  fact  T^  and  T^  are  not 
associated.  This  follows  from  the  fact  (Johnson  and  Kotz  (1972),  p.  262)  that 
they  are  negatively  correlated. 

4. 2.  The  multivariate  case. 

Let  ]_  =  (T^,...,T  ),  n  >  3,  be  a  nonnegative  absolutely  continuous 
random  vector  as  in  Section  3  and  let  T_  be  defined  as  in  (3.3)  and  (3.4). 

As  in  Section  4.1  we  will  not  distinguish  between  T  and  T  and  just  write 
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To  clarify  the  general  condition  in  Theorem  4.4  below,  consider  the 
case  n  =  3.  Following  (3.5)  -  (3.7),  for  x.  =  (x^.x^jX^)  which  satisfy 
aJ^Xj)  <  min  (a^1  (xg  ),h^1  (x3))  and  A^O^A^A^1^))!  Tj  ■  Aj^Xj)) 

^  A3  (X3-A3 (a^  (x^))|Tj  -  A^  (x,))  let 

ai (X1.X2.X3)  =  A ^  (x^), 

a2  (X1  ,x2,x3  )  =  a  1  00  A2  [  x2  ”  ’I' 2 1  1  (al  OS  ^  I  a  1  00 )  |  T 1  _  a ^  (xj ] 

=  ax(x)  +  A2  [x2  -  A2(a1(x))|T1  =  aj(x)], 

a3  ^X1 ,x2 ’x3  ^  =  a2 00  +  ^3  [ x3"^3| 1  2^a2^0l  a l  ^ ’ a 2 ^ ) I  ^ 1  ~  al 00 5 ^2  ~  a2 ^  1 

*  _  j 

=  a 2 00  +  A3  [x3  -  A3(a1(x) ) 

_  A 3C a 2^ii)  -  —  a^(2<))|T^  _  a^(^0,  T2  - 


Clearly,  on  the  given  domain,  each  ai  increases  in  xj,  each  of  83  and 
33  increases  in  X2  and  83  increases  in  x3.  Thus,  to  find  conditions  for 
association  of  Tj,  T2  and  T3  we  only  need  to  find  conditions  such  that 

(4.9)  a2  increases  in  x^, 

(4.10)  83  increases  in  X2> 

(4.11)  a3  increases  in  xj. 

As  in  Section  4.1  it  can  be  shown  that  (4.9)  is  the  same  as 
^2)  1  (t  2!  1 1 )  >  >i'2j  1  (tjl 1  i+u^’  t2>tl  +  u>tl  ^  or’  POre  exPlicit1y> 

(4.12)  A2  (t  ^  )  +  A2  (t2  “t  1 1  T1  =t  ^  )  >  A2(t1+u)  +  A?  ( 1 2  ~t  1  -u|  Tpt^u),  t?  >  tj+u  >  tj  1 


or,  by  means  of  derivatives. 


X2(t1)  c  x2  ^  t2 1  T1  tl’*^  ”  acj  A2(t2"tl!TrCl1|c1=t1’  l2  *  tl 


0.  Similarly, 
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(4.10)  is  the  same  as 


*3|  1  2^3!  ti,t;2^  >  ^3)1  2(t3|  t1,t2+u),  t3  >  t2  +  u  >  t2  >  t1  >  0,  or  more 
explicitly, 


(4.13)  ^3 ( 1 2 ~t  1 1  T l ~t ^  )  +  A3  ( t ^  —  1 2 j  T^-t^  *  "^2  =  ^2^  * 

^3 ( 1 2  +u ~t 1 1 T 1 =t 1  )  +  ( 1 2 ~t2 _u | T| -t ^ ,  T2-t2+u ) ,  t3  >  t2  +  u  >  t2  >  t^  >  0, 

or,  by  means  of  derivatives, 

^  3  ( 1 2 1  T 1 =t  1  *  * )  <  x3(t3lTi=tpT2=t2,*)  ”  !cj  A3^t3't2lTl=tl’T2=C2^C2=t2  ’ 
t3  >  t^  >  t^  >  0  .  Finally,  (4.11)  is  the  same  as  (here 

^2  (*2,a  )  h  a  +  A^1 (X2"A2 (a )|  T1 =a ) , x2  >  °»  a  >  °) 

4* 3 j  1  2 C 1 3I  tpa2(x2*t  j))  >  3 J  1  2^31^1"*"*^ > a  2  ^  *  2  ’ ^  1  ^  ’ 

x2  >  0,  0  <  t1  c  t1  +  u  <  a2(x2,t1)  <  a2(x2,t1+u)  <  t3  ,  or,  more  explicitly, 

Aj(t  1)  +  A3(a2(x2»t  p  -  t1jT1  -  t^)  +  A3(t  3~a  2(^2*t  ^  )|  }_t  2(^2»  t  ^ ) ) 

(4.14)  ?  A3(t1+u)  +  A3(a2  (x2,t^+u)-t1-u|  T^t^+u) 

+  A3(t3-a2(x2.t1+u)|T  j-t^u,  T 2 _a  2 ( x 2 » t  ) )  , 

^2  ^  0 ,  0  ^  1 ^  <  t  ^  u  ^  ^2^2* ^1^  ^  a3(x2,t  ^+u )  ^  t3  , 

(see  Figure  4.2)  or,  by  means  of  derivatives. 


«  •.  ■  •*.  /  ,* 


•  '  *  .  *  .  ■  .  .  .  *  ,  .  .  «  ■  •  •  ■. 

.  .  *  *  *  *  V-  *  -  *  .  1  •  V  •  '-f  ‘ 


",3|l,2(t'|tl-C2) 

!  /  =A3(ti^+A3(c2-tl^Tl=tl^ 

|/  +A3(t-t2|T1=t1,T2=t2) 


lJj3|l,2(t,tl+U*E2) 

=A3(t1+u) 

/  +A3(t2-t1-u|T1=t1+u) 
+A3(t-t2|T1=t1+u,T2=t2) 

^3 | 1 (t 1 C1^=A3 (tl^ 

>  +A3(t-t1|T]=t1) 


’^3l  l(t!tl+u)=A3(tl+u) 
+A3 | T^t^+u) 


A  3  ( t ) 


t2=32(x2’tl) 


1 2=a9 (x2 , t 3+u) 


Figure  4.2.  Illustration  of  (4.14). 
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(4.15) 


3a?(x?,a) 

x3 (t  1 )  +  x3^t2lTl=tr*  ^  3a 


+  3c,  A3(t2“tllTl“cl>l  _t 

a=ti  1  j  c  i  =t  i 


+  ac,  A3(t3_t2lTl=cl*  T't2)| 


cl=ti 


3a2(x2,a)| 


+  ac?  a3 (t3_t2l Tl=tl *  T2=c2 M  *  3a  . 

^  I c2~z2  j 

3a?(x?,a) 

*  X3(t2|T1-t1,. )  +  x3(t3|T1=tpT2=t2)  .  — — 


where  in  (4.15)  t2  =  a2(x2,t^). 

In  the  statement  of  the  next  theorem  we  use  the  notation  (see  (3.4)) 

*k  (xk  *ti  •  •  •*»tk-l  ^-tk-l+Ak  f  xk‘^k|  1, . . . ,  k-1  ^k-ll  tl . tk-l  ^1  T1  =tl  ’  *  *  *  ,Tk-l  =tk- 

which  describes,  according  to  (3.4),  how  is  determined,  given  that 

A  A 

j  ^  j|(  —  ^  and  t  h  3 1  X  ^  —  *  T  ^  _  ^  —  t  ^  ^  • 


Theorem  4.4^_  If  1>2 . . tk. 


^k*l^xk+l,tl,"-’tk*’ak,2^xk*2,tl . tk,“k+l* . 


'  k ’  k+1 ‘ 


®i-l(xi_l»t]i».*.,ti<,  a|c+i,...,a’i_2)]  (the  arguments  of  some  of  the  a^’s 


are  omitted)  decreases  in 


tk  e  (tk * t k  ^  tk-l*  ak+l  ^  ^k’  ^t+1  ^  ^  k+1, . . . , i -lj  for  all 

1  <  k  <  i  <  n,  x^  >  0,...,x.  >  0,  0  <  tj  <  t2«...<tk_p  and  if  the  above 
condition  holds  for  all  permutations  of  the  indices  1,2, ...,n,  then 


Tj,...,Tn  are  associated. 

Sketch  of  the  proof.  The  conceptualy  simple  but  notationally  involved  proof 
of  Theorem  4.4  is  similar  to  the  proof  of  Theorem  4.1  and  we  will  omit  the 
details.  The  idea  of  the  proof  is  as  follows:  According  to  (3.4),  given 
that  ~  l,*«»,j-j_^  ~  i  —  1,  j-j  “  i  and  that  T^  —  t^,...,T,|  ^  —  t  ^  ^ 

=  x.j ,  the  value  of  T^  then  is 


and 


Equivalently,  is  the  solution  of 


(4.16) 


x.  =  tl>. 

l 


i) l,...,i-l'  i‘ 


|  ti » • • • ) 


Association  of  Tj,...,Tn  will  follow  if  we  show  that  t^  increases  in  xk(k 
<  i).  In  order  to  do  that,  we  can  fix  ti,...,tk_i  and  express  tk  as  a 
function  of  xk  and  express  tk+}, • * •>ti _i  as  functions  of  tk  (and  of 
xk+1,...,xi_1).  Since,  trivially,  tk  increases  in  xk,  it  suffices  to  show 
that  ti  increases  in  tk  (recall  that  t^,...,tk_j  and  xk,...,x^_1  are 
held  fixed). 

Rewrite  (4.16)  as 


(4.17) 


xi 


k-i 

W  ‘  L  Ai (Wil Ti=ti . VrVi^ 

1-2 


£  Ai(VVilTrti,*,,’Vrt£-i) 

=K 

+  Ai (ti-ti_1| T1=t1,...,Ti_1=ti_1) 


and  notice  that  t7-  is  determined  as  the  solution  of  (4.17).  Clearly 


(4.18)  (u| T^=tj,.. .,"T  =t^_| )  increases  in  u  >  0  . 

The  LHS  (4.17)  is  fixed  (i.e.,  it  does  not  depend  on  tk).  If  the  condition 
of  Theorem  4.4  holds,  that  is,  if  for  a  fixed  t^,  RHS  (4.17)  decreases  as 
tk  increases,  then  (using  (4.18))  the  solution  ti  of  (4.17)  must  increase 
as  tk  increases.  But  this  is  what  we  wanted  to  prove,  >i 


Remark  4.5.  Theorem  4.4  should  be  contrasted  with  Theorem  4.7,  p.  146,  of 
Barlow  and  Proschan  (1975).  They  show  that  if  for  all  t^. 


(4.19)  P(T.  >  ti|T1 


ti , . . . ,Ti _i 


ti-l) 


increases  in  t 


r 


*  *  * ti -1  * 


i  =  2, 


then  -  Tj, .. . ,Tn  are  associated.  Their  proof  essentially  constructs 

Tj,...,V  as  in  (3.10)  and  (3.11)  and  then  argues  that  (4.19)  implies  that 

T',...,V  are  increasing  functions  of  Uj,...,Un  of  (3.10)  and  (3.11).  In 

proving  our  Theorem  4.4  we  follow  the  same  line  of  thought  but  apply  it  to 

T, . T  which  arise  from  the  total  hazard  construction  described  in  (3.3) 

1  n 

and  (3.4). 


Remark  4.6.  Shaked  and  Shanthikumar  (1984b)  showed  that  if  for  disjoint  sets 
I,  J  C  {l,...,n}  and  fixed  t_j,  _tj,  tj  [such  that  t_j  <  t  j  ]  and  k  e  TOU 
(I  or  J  may  be  empty), 

(4.20)  XkC(  V  t.)  V  (  V  t.)  +  u|Tj  =  tj,  Tj  =  tj,.) 
iel  *  jeJ  J 

>  Ak((  V  t.  )  V  (  V  t.)+u|Tj  =  tj,.),  u  >  0, 
iel  jeJ  J 

then  T^, . . . ,Tn  are  associated.  Below  it  is  argued  that  (4.20)  implies  the 
conditions  of  Theorem  4.4.  Thus  Theorem  4.4  provides  a  new  route  of  proving 
Theorem  5.2  (and  Remark  5.5)  of  Shaked  and  Shanthikumar  (1984b). 

To  avoid  messy  notation  we  consider  the  case  n  =  3  and  show  that  (4.20) 
implies  (4.12),  (4.13)  and  (4.14). 

Proof  that  (4.20)  ==>  (4.12).  From  (4.20)  we  get 


26 


(4.21)  x2(t1+v|.)  <  x2(t1+v|  Tj=tj,» ),  t-^  >  0,  v  >  0. 


Integrating  (4.21)  with  respect  to  v  over  [0,u]  we  obtain 


(4.22)  A2(t^+u)  -  A2 (t ^ )  <  A2 (u|  T^=t^ ),  u  >  0,  tj  >  0. 


From  (4.20)  we  also  get 


(4.23)  x2  (t^+u+vj  T^=t^+u,  •  )  <  X2(t^+u+vj  Tj=tp« )»  t^  >  0,  u  >  0,  v  >  0. 


Integrating  (4.23)  with  respect  to  v  over  [O^-tj-u]  we  obtain 


(4.24)  A?  (t2-t^  -u|  T  j  =t  j+u )  <  A2  (t2-t  y\ T1=t1 )  -  A2(u|T1=t1),  t2  >  t^u  >  t1  >  0. 


Adding  (4.22)  and  (4.24)  and  rearranging  one  obtains  (4.12). 


Proof  that  (4.20)  ==>  (4.13).  From  (4.20)  we  get 


(4.25)  *3^2+Vl  ^l=tl’* )  <  x3^2+vl  ^l=tl’^2=t2’*  *1  >  0*  *2  *  u  >  v  >  0* 


Integrate  (4.25)  with  respect  to  v  over  [0,u]  to  obtain 


(4.26)  A^  (t2+u-t  J  T^=t ^  )  "  43(^2 -ti!  Ti=ti  ^  4  Ag (u|  =t ^ » T2  =t^  ), t ^  +u> t2>t|>0. 


Condition  (4.20)  also  yields 


(4.27)  X3(t2+u+v|  T1=t1,T2=t2+u,*  )  <  X3  (t2+u+v|  =t^  ,T2=t2,  •  ),  t^>  0,  t2>  0,  u>0,  v>0. 


*  •  .'*•*-*  - v  »V.‘»  *’  .  *  .  *  .  *  . 


Integrate  (4.27)  with  respect  to  v  over  [0, 1 3-t 2-u 3  to  obtain 


(4.28)  A3 (t 3 -t2-u j T1=t1,T2=t2+u)  <  A3 (t 3 -t2|  Ti=ti»T2=t2) 

-  A3(u|T1=t1,T2=t2),  t3  >  t2  +  u  >  t2  >  tj  >  0. 

Add  (4.26)  and  (4.28)  and  rearrange  to  obtain  (4.13).  11 

Proof  that  (4.20)  ==>  (4.14).  The  following  follow  from  (4.20)  [See  Figure 
4.2]: 

(4.29)  X3(t1+v|.)  <  X3(t1+v|T1=t1,.),  tj  >  0,  v  >  0, 

(4.30)  X3(t^+u+v|  Ti=tj+u,» )  <  X3(t1+u+v|  T^tp* ),  t^  >  0,  u  >  0,  v  >  0, 

(4.31)  X3(a2(x2*t  ^)  +  v |  T  jat  j+u , • )  <  x3(a2(x2,t1)+v|T1-t1,T2-a2(x2,t1),. ), 

t^  >  0,  u  >  0,  v  >  0,  x2  >  0, 

(4.32)  x3(a2(x2,t1+u)+v| T^j+u,  T2=a2(x2,t1+u),«) 

^  X3  (a2  (x2,tj  +u )  "^  v  j  T  ^  —  t  j,T2-a2(x2,tj  ),  • ),  t  ^  >  0,  x2  >  0,  u>  0,  v  >  0. 

Integrate  (4.29)  with  respect  to  v  over  [0,ul  to  obtain 

(4.33)  A3(t1+u)  -  A3(t1)  <  A3(u|T1=t1),  0  <  tj  <  tj  +  u. 

Integrate  (4.30)  with  respect  to  v  over  [0,  a2(x2,t^)  -  t^  -  u]  to  obtain 

(4.34)  A 3 (a 2 ( x 2, t ^ ) ~t | ~u j T ^ =t ^+u )  <  A3(a 2( x2,t ^ ) -t T ^ -t 

*  A 2 ( ^ |  ^ — T ^ ) *  0  ^  t  ^  ^  t  ^  u  <  a2(x2»tj). 

Integrate  (4.31)  with  respect  to  v  over  [0,  ^(x^.t  i+u )-a ?( x^.t  ^  ]  to 


obtain 


(4.35)  A3(a2(x2,t1+u)-t1-u(T1=t1+u)  -  A3(a 2(x2,t  1  )-t ^uj T^t ^u) 

*  A3  ( a2  ( x2  *  ^  i +u  )  "a2  ( x2  •  1 1  )|  T  j  =t  j ,  T2  =a£  ( X2  »t i ) ) , 

0  <  t2  <  t2  +  u  <  a2(x2,t2)  <  ^(^jti+u). 

Finally  integrate  (4.32)  with  respect  to  v  over  [0,  t3  -  a2 (x^t^+u) ]  to 
obtain 

(4.36)  A3(t3-a2(x2»t1+u)|T1=t1,T2=a2(x2,t1+u)) 

<  A3(t3-a2(x2,ti)|Ti=ti »T2=a2  (x2  *  ^ l  ^ ) 

-  A3(a2(x2,t  j+u )-a2(x2»t^)|  Tj-t^,T2-a2(x2,t  ^ )) , 

0  <  tj  <  tj+u  <  2*2 (x2 * 1 2 )  <  a2 ( x2 » 1 2  +u )  <  t3. 


Add  (4.33)  -  (4.36)  and  rearrange  to  obtain  (4.14).  11 


5.  An  application:  stochastic  ordering. 

5.1.  The  bivariate  case. 

Let  (^  1 » S 2 )  and  (Ti^)  be  two  nonnegative  absolutely  continuous 
random  vectors.  The  corresponding  hazard  rates  and  cumulative  hazard 
functions  will  be  denoted  as  follows: 


(t) 
ri  (t) 
Q1  (t ) 

Ri  (t) 


1 im  — p  P{t  <  S,  <  t  +  At |  S, 
AtiO 

lim  — r-  P(t  <  T.  <  t  +  At j  T , 
AtiO 

/ 0  Mi  (u)du,  t  >  0,  i  =  1,2, 
/q  r i  (u  )du ,  t  >  0,  i  =  1,  2, 


>  t ,  s2 

>  t ,  T2 


q x (t ) 1 2 )  =  lim  ~  P { t  <  S1  <  t  +  At | S ^  >  t. 


>  t),  t  >  0,  i  «  1,2, 

>  t),  t  >  0,  i  =  1,2, 

^2  =  ^  t  5  ^2  > 


1 


q2(t|tx)  = 

•”  1  (^ )  ^2  ^  = 
r2(tlt1)  = 


=  lim  —r  P{t  <  S2  <  t  +  At | Sjat ,,  s2  >  t},  t  >  t,  >  0, 
At-vO  * 

=  lim  P{t  <  Tj  <  t  +  At ) Tj  >  t,  T2  =  t2) ,  t  >  t2  >  0 
At-t-0  , 

=  lim  -r  P{t  <  Ty  <  t  +  At|Ti=t,,  To  >  t} ,  t  >  t-,  >  0, 
At+O 


Q1(t| t2 )  = 


o2(t|t1)  = 


R1(t|t2)  = 


R2(t|t1)  = 


t?+t 

Jt  (u | t2)du,  t  >  0,  t2  >  0, 
t,+t 

lt |  q2(uj  tj  )du,  t  >  0,  tj  >  0, 
t2+t 

/t  r^(ujt2)du,  t  >  0,  t2  >  0, 
t.+t 

ft  r2 (u 1 1 1  )du,  t  >  0,  tj>  0. 


The  inverses  of  Q.(.),  Q  -j  ( *  1 1 ) ,  R^O  and  R  i  ( •  1 1 )  are  defined  in  an 
obvious  manner  as  in  Section  2. 

Let  Xj  and  X2  be  independent  standard  exponential  random  variables, 
Define  Sp  S2  as  follows:  On  {O^(X^)  <  021(X2)}  let 


(5.1)  S:  =  Qj1 (Xj), 

(5.2)  S2  =  Q j1 (X x )  +  Q”1CX2-Q2(Qi1 (xi) )j  0"1 (Xx )), 


and  on  (Q^(X^)  >  ^2^ (X2 ^ 


(5.3)  Sj  =  Q^1  (X2)  +  Qj1(XrQi(Q21(X2))|Q21(X2)), 

(5.4)  S2  h  021(X2). 


Similarly  on  {Rj^(X^)  <  R2^(X2)}  let 


*•  *•  *.  ■'w  ~  “  *  '  * . '  •  •  "  »  *  -  ‘  •  ' 

J1  *  -  .*  •  -  V  •  *  *  A  - 


xmpuu  mm aaetaat, 


a 


(5.5)  Tj  =  RJ1(X1), 


(5.6)  T2  =  R~1(X1)  +  R21(X2-R2(Rj1(X1))|Rj1(X1))> 


and  on  {Rj^(Xj)  >  R^1^)}  let 


(5.7)  Tx  =  R^1  (X2)  +  Rj1(XrRi(R21(X2))|R21(X2)), 

(5.8)  T2  e  R2!(X2). 


Theorem  5.1.  Let  (S^,S2)  and  (T^,T2)  be  defined  as  in  (5.1)  -  (5.8) 
where  X^  and  X2  are  independent  exponential  random  varibles.  Then 

(slfs2)  (slts2), 

(T1#f2)  &  (JVJ2). 

Proof:  Apply  Theorem  2.1  twice.  H 

In  (5.1)  -  (5.8)  we  use  the  total  hazard  construction  twice:  to  construct 

A  A 

X  and  to  construct  X  •  Note  that  we  use  the  same  X^  and  X2  Tor  both 
constructions.  Thus,  roughly  speaking,  we  put  X  and  X  on  the  same 
probability  space.  This  enables  us  to  compare  them  realization-wise  as  is 
done,  e.g.,  in  the  next  theorem. 

Theorem  5.2.  If 


(5.9) 

(5.10) 

(5.11) 

(5.12) 


Qj(w)  >  Rj (w),  w  >  0, 

02(w)  >  R2(w),  w  >  0, 

*^2  l )  +  ^2  ^w_s  1 1  s  1 )  >  R2(w).  w  >  s1  >  0, 

Ql (s2 )  +  Q1  (w-s2|s2)  >  Rj(w),  w  >  s2  >  0, 


%  •  *  »  1  C"  -  ■«  “  »  ’  •  *»  «  ■  »  .■  •_*  -  x  *  I  '  -  •  •  "  *  '  •  *•“„*.*, 

•j  V  '.*  *  V  '  -  V  V  V  V  v  V  j  *.•  *  •  v>  '  -  V  S  1  . 


V  •  .*•  N  **  .*•  .% 

•  .*  W* 
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(5.13)  (^(Sj)  +  02  (w~si  J  s  i )  >  R2  (t  x  )  +  R^w-tjjtj),  w  >  tj  >  Sj  >  0, 

(5.14)  Qj^)  +  01(w-s2|s2)  >  (t2 )  +  R1(w-t2|t2),  w  >  t2  >  s2  >  0, 

then 

(5.15)  (SltS2)  f1  (Tj.Tg). 

Remark  5.3.  The  conditions  of  Theorem  5.2  simply  state  that  at  any  time  w 
(no  matter  what  the  previous  history  is)  the  cumulative  hazard  of  Sj  is 
larger  then  the  cumulative  hazard  of  T^,  i  =  1,2.  The  proof  of  Theorem  5.2 

A  A 

below  uses  the  fact  that,  since  the  total  cumulative  hazards  of  S .  and  T . 
by  the  time  they  failed  must  be  equal  (to  X^),  then  necessarily  <;  7\ 
a.s.,  i  =  1,2  (see  Figure  5.1  for  a  typical  realization  when 

Q^lXj)  <  q21(x2)). 

Proof  of  Theorem  5.2.  We  will  show  that 

A  A 

(5.16)  <  T^  a.s.  , 

(5.17)  S2  <  T2  a.s.  , 

and  the  result  then  follows  from  Theorem  5.1  and  (1.5). 

Let  X}  and  X2  be  independent  standard  exponential  random  variables. 
First  consider  the  case  Qj^(Xj)  <  02^(X2).  Then 

sx  =  oJ1(x1)  , 

(5.18)  S2  =  Sj  +  Q21(X2‘Q2(^l)iS\)  * 


If  RJ^Xj)  <  then  fj  =  Rj1^)  >  Q*1^)  =  Sj  and  (5.16) 

(5.9) 

follows.  Also,  in  this  case 

(5.19)  f2  =  f,  +  R-1(X2-R2(f1)|f1)  . 

Thus-  T9  =  inf{w  >  f1:R9(T,)  +  R9(w-T1|f1)  >  X9) 

2(5.19) 

>  inf{w  >  S.  :Q0 (S, )  +  Q0(w-S1|S1)  >  X9}  =  S,  and  (5.17)  follows. 

(5.13)  121  2  1  1  2  (5.18)2 

If  Rj^Xj)  >  R2V2)  then 

(5.20)  T2=R21(X2), 

(5.21)  fj  =  f2  +  R^1(X1-R1(T2)|f2)  . 

A  A 

By  assumption,  Sj  <  S2  .  Thus, 

(5.22)  Sj  <  S2  =  inf{w  >  SjiQgfSj)  +  Q^w-Sj^)  >  X?} 

(5.18) 

<  inf{w  >  S,  :R9(w)  >  X9}  . 

(5.11)  1  2  2 

A  A  1 

Also,  R2^S1 )  (5<10)(^2^S1  ^  =  ^2^1  ^Xl^  <  X2  ’  where  the  second  inequality 
follows  from  the  assumption  Q^(X^)  <  ^2  ^  ( x  2 )  •  ^hus*  since  R2(w)  increases 
in  w  >  0, 

(5.23)  inf (w  >  Sj :R2 (w)  >  X2)  =  inf{w  >  0:R?(w)  >  X2) 

A  A 

^2  <  ^ 1  * 

(5.20)  (5.21) 1 

A  A  A  A 

Combining  (5.22)  and  (5.23)  it  follows  that  <  S2  <  T2  <  T^  and  both 
(5.16)  and  (5.17)  follow. 

The  proof  of  (5.16)  and  (5.17)  for  the  case  Q^X^)  >  0^ 1  ( x 2 )  is 


similar,  using  (5.10),  (5.19),  (5.12)  and  (5.9)  instead  of  (5.9),  (5.13), 
(5.11)  and  (5.10)  respectively,  n 

5.2.  The  multivariate  case. 

Let  =  (S^,...,Sn)  and  T  =  (T^,...,Tn)  be  nonnegative  absolutely 
continuous  random  vectors.  Denote  by  q.  (s|  Sj  =_s^,*)  the  conditional  hazard 
rates  of  defined  as  in  (3.1)  and  denote  by  Q ^  ( s | =  s^  )  the  conditional 
cumulative  hazards  of  defined  as  in  (3.2).  Similarly  denote  the 
conditional  hazard  rates  and  cumulative  hazards  of  T_  by 

ri  (tllj  =  and  Ri^llj  = 

Using  the  total  hazard  construction  (see  Section  3)  one  can  express 
S  and  J  (such  that  S  ^  S  and  T  T_)  as  functions  of  the  same 
independent  standard  exponential  random  variables  Xj,...,Xn.  Using  these 

A  A 

S_  and  _T  one  can  prove  the  following  result  using  the  method  of  the  proof 
of  Theorem  5.2  (but  with  more  involved  notation).  We  omit  the  details. 

Theorem  5.4.  If  for  1  <  i  <  j  <  n,  0  <  s^  0  <  t^  «...•;  t^  , 

0  <  s^  <  t..,  i  =  1,...,£,  and  all  permutations  tt_  of  (l,...,n), 

(5.24)  Qk(si+rSilSn(l)  =  Sl’**',Sn(i  )  =  Si  ) 

£_lk  (w_sj!  Sif  (1)  =  S1 . Sir(j  )  =  Sj) 

>  Rk  ^i+l _ti I  Ttt ( 1 )  =  tl»  —  ,Ttt (i )  =  ti  ) 

+  Rk(w-t)l|TiT(1)  =  ^1,---,Tir(£)  *  V 


whenever  w  >  s.  V  t  (where  empty  sums  are  identically  zero)  then 


3' 


The  proof  of  Theorem 

5.4  consists  of  constructing 

S_ 

and 

T 

by 

(3.3) 

(3.4)  using  the  same 

X j, . . . ,Xn  (that  is,  putting 

S 

and 

T 

on 

the  same 

probability  space)  and  noticing  that  (5.24)  implies  that  realization-wise 
_S  <  T  .  The  result  then  follows  from  (1.5). 

Remark  5.5.  Using  the  standard  construction  (3.10)  and  (3.11)  one  can  show 
the  following  analog  of  Theorem  5.4  (see,  e.g.  Veinott  (1965)  or  Arjas  and 
Lehtonen  (1978)):  If 

(5.26)  Sj  Tj_ 

and  for  <  t^,...,s^_^  <  t._j, 

(5.27)  [S-jl^i  -  Sj,...,S^_j  -  1  [TjTi  =  t|,...,T^_^  = 

i  =  2,3,...,n,  then  _T.  The  idea  of  the  proof  of  (5.26)  +  (5.27)  ==> 

(5.25)  is  the  sane  as  the  proof  of  Theorem  5.4:  Using  the  same  U^,...,un  of 
(3.10)  and  (3.11),  put  _S  and  T_  on  the  same  probability  space  and  note  that 

(5.26)  and  (5.27)  imply  that  real i zati on-wi se  S_'  <  T'  . 

Remark  5.6.  Shaked  and  Shanthikumar  (1934b)  proved  that  if  for  all  disjoint 
sets  I,  J  C  {l,...,n}  such  that  Tim  t  |d  and  for  all  fixed  v,  ;>  0e  the 
following  holds: 


Zi*  X)’’ 

>  '•kMl|  -  u  >  l.v  »j)  v  <_v  »(>• 

JeJ  1 e I 


A  \W  \* 


(5.28) 
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whenever  >  _v j »  u  >  0  and  k  e  I  U  J  (I  or  J  may  he  empty)  then 

c  st  T 
S  <  T  . 

Theorem  5.4  provides  a  new  way  of  proving  their  result.  In  fact  Theorem 

5.4  is  a  stronger  result  than  Theorem  3.1  and  3.4  of  Shaked  and  Shanthikumar 

(1984b)  because,  as  will  be  argued  shortly,  (5.28)  implies  (5.24). 

•To  see  that  (5.28)  implies  (5.24)  suppose  that  in  (5.24),  jt_  =  (l,...,n) 

and  order  the  j  s/s  and  z  t^'s  of  (5.24)  in  an  increasing  order 

0  <  v,  <...<v.,  where  each  v_  is  either  s,  or  t.-  for  some  i.  Also 

1  j+Jt  f11  1  1 

set  vQ  =  0  and  Vj+ji+1  =  w  (the  w  of  (5.24)).  Consider  some  interval 
EVvm+l)*  Let  1  =  {i:W*  J  =  *iEl’  si  <  V  •  Note  that  if 

i  e  I  then  s^  <  v^  because  in  (5.24),  s^  <  t^  for  i  =  1,...,£  .  Let 

u  e  tvm,vm+1)  and  apply  (5.28)  with  the  above  I  and  J  and  _Vj  =  s^  , 

A 

*!  =  tj,  Vj  =  Sj  to  obtain 

(5.29)  qk^ul^-I  >  rk  (u)  Xi  =  li»-)  • 

Integrating  (5.29)  with  respect  to  u  over  [v^.v^+j)  and  adding  the 
resulting  integral  inequalities  over  m  =  0,...,j  +  z  one  obtains  (5.24)  with 
jt_  =  (1,2, ...,n).  The  proof  for  other  permutations  £  is  similar. 

For  example,  to  obtain  (5.9)  note  that  (5.28)  implies 

(5.30)  q^(uj.)  <  r^(u)*)»  u  >  0. 

Integrate  (5.30)  with  respect  to  u  over  [0,w)  to  obtain  (5.9).  To  obtain 
(5.11)  note  that  (5.28)  implies 

(5.31)  q^ujSj  =  Sj,*)  <  (uj  • )»  u  >  Sj  . 


Integrate  (5.30)  with  respect  to  u  over  [0,s^)  and  integrate  (5.31)  with 
respect  to  u  over  [s^,w)  and  add  the  two  resulting  integral  inequalities  to 
obtain  (5.11). 

6.  Further  applications. 

6.1.  tiulti variate  increasing  failure  rate  average  (MIFRA)  distributions. 
Theorems  4.1,  4.3  and  4.4  give  conditions  under  which  the  distribution 

of  Ti,...,Tn  can  be  expressed  as  the  distribution  of  increasing  functions  of 
independent  exponential  random  variables.  If  these  increasing  functions  are 
also  subhomogeneous  (a  function  g:  R^  -*■  R+  is  subhomogeneous  if 
g(at)  <  a  gU)  for  all  a  e  [0,1],  t  >  0,  see,  e.g.,  Marshall  and  Shaked 
(1982))  then  _T  satisfies  the  MIFRA  condition  of  Block  and  Savits  (1980). 

For  example  in  (2.21)  [when  a  <  a',  B  <  B‘]  and  are  expressed  as 

incresing  subhomogeneous  functions  of  Xj  and  X£.  Hence  (T^^)  of  Example 
1  (when  a  <  a',  B  <  B1)  is  MIFRA.  This  result  has  been  obtained  also  by 
Marshall  and  Shaked  (1982)  and  Shaked  (1984).  A  special  case  of  this  result 
can  be  found  in  Block  and  Savits  (1980). 

6.2.  Variance  reduction  in  simulation  of  dependent  variables. 

Let  S.  =  (Sj,...,S  )  and  T  =  (Tj,...,Tn)  be  random  vectors  and  let 
g:Rn  -*■  R  and  h:Rn  R  be  monotone  in  the  same  (or  the  opposite) 
direction.  Due  to  theoretical  or  technical  reasons,  the  expected  value 

(6.1)  E[g(S)  -  h(T)] 


may  be  hard  to  compute  in  some  applications.  One  possible  recourse  is  a 
simulation  of  g(S)  and  h(_T).  That  is,  k  independent  replications  of 


and  ]_  are  generated  using  pseudo-random  numbers,  and  (6.1)  is  then 
estimated  by  averaging  the  k  realizations  of  g(Sj  -  h(T). 

Rubinstein,  Samorodnitski  and  Shaked  (1985)  have  considered  an  efficient 
method  of  simulating  g(Sj  and  h(]T)  when  the  distributions  of  and  T 
satisfy  (4.19).  Their  method  is  based  on  the  fact  that  when  (4.19)  holds, 
then’  g(Sj  and  h(^T)  can  be  represented  as  increasing  functions  of 
independent  uniform  [0,1]  random  variables,  using  the  standard  construction 
(3.10)  and  (3.11).  Then,  by  putting  and  T_  on  the  same  probability 
space,  one  can  reduce  the  variance  of  the  Monte  Carlo  estimate  of 
E(g(S)-h(T)). 

The  same  idea  may  apply  for  random  vectors  with  distributions  satisfying 
the  condition  of  Theorem  4.4.  Under  this  condition  too  it  is  possible  to 
represent  g(Sj  and  h(T)  as  increasing  functions  of  independent  random 
variables,  put  them  on  the  same  probability  space  and  reduce  the  variance  of 
the  Monte  Carlo  estimate. 


In  some  applications,  even  if  both  (4.19)  and  the  condition  of  Theorem 
4.4  hold,  the  total  hazard  construction  (3.3)  and  (3.4)  may  yield  simpler 
expressions  than  the  standard  construction  (3.10)  and  (3.11).  In  such  cases 
use  of  the  total  hazard  construction  is  prefarable.  A  study  of  these 
ramifications  of  the  total  hazard  construction  is  planned. 

6. 3.  Multi -unit  imperfect  repair. 

Shaked  and  Shanthikumar  (1984a, b)  considered  a  model  for  imperfect  repair 
of  multi-unit  systems.  In  that  model,  n  units  (whose  original  lives 
T^,...,Tn  have  absolutely  continuous  distribution)  start  to  live  at  the  same 
time.  Upon  failure  an  item  undergoes  a  repair  and  is  scrapped  if  the  repair 
is  unsuccessful.  If  i  items  (i  =  0,l,...,n-l)  have  already  been  scrapped. 


then,  with  probability  pj+j  the  repair  is  successful  and  the  item  continues 
to  function  -  but  it  is  only  as  good  as  it  was  just  before  the  repair  -  and 
the  other  items  "do  not  know"  about  these  failure  and  repair.  With 
probability  1-P-j+i  the  repair  is  unsuccessful  and  the  item  is  scrapped. 

Mathematically,  if  the  original  lives  have  the  conditional  hazard  rates 
AkHli  k  e  T,  then  the  resulting  lives  Tlt...,T  have  the 

conditional  hazard  rates  xk  given  by 

(6.2)  AI((*|^I  =  .tj » • )  =  P|  i  j  +1  A|<  ( ‘  I  Xj  ^  e  I* 

where  |l|  is  the  cardinality  of  I. 

From  (6.2)  it  follows  (see  (3.2))  that  Ak(t|Tj  =  _tj),  k  eT  —  the 

hazard  accumulated  by  Tk  during  the  time  interval 

(V  t.,  V  t .  +  t],  t  >  0,  --  is  given  by 
i  el  1  i  el  1 

(6.3)  =  j )  =  P|l|+l  \  ^1  —I  =  —I  ^  e  ^  • 

Using  (6.3),  various  results  of  this  paper  can  be  restated  for 
Ti,...,Tn.  For  example,  if  p^  =  pg  =...=pn  and  the  a^'s  satisfy  the 
conditions  of  Theorem  4.4  then  Tp...,Tn  are  associated.  Similarly,  two 
random  vectors  resulting  from  application  of  imperfect  repair  can  be  compared 
stochastically  if  the  original  random  vectors  satisfy  the  conditions  of 
Theorem  5.4.  Proposition  6.6  of  Shaked  and  Shanthikumar  (1984b)  can  be  proven 
using  these  ideas. 
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